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ABSTRACT 

Energy release in the early Universe (z < 2 x 10 6 ) should lead to some broad spectral distortion of the cosmic microwave background 
(CMB) radiation field, which can be characterized as v-type distortion when the injection process started at redshifts z < 5 x 10 4 . Here 
we demonstrate that if energy was released before the beginning of cosmological hydrogen recombination (z ~ 1400), closed loops of 
bound-bound and free-bound transitions in H i and He n lead to the appearance of (i) characteristic multiple narrow spectral features 
at dm and cm wavelengths, and (ii) a prominent sub-millimeter feature consisting of absorption and emission parts in the far Wien 
tail of CMB spectrum. The additional spectral features are generated in the pre-recombinational epoch of Hi (z > 1800) and Hen 
(z ^ 7000), and therefore differ from those arising due to normal cosmological recombination in the undisturbed CMB blackbody 
radiation field. We present the results of numerical computations including 25 atomic shells for both H i and He n, and discuss the 
contributions of several individual transitions in detail. As examples, we consider the case of instantaneous energy release (e.g. due to 
phase transitions) and exponential energy release because of long-lived decaying particles. Our computations show that due to possible 
pre-recombinational atomic transitions the variability of the CMB spectral distortion increases when comparing with the distortions 
arising in the normal recombination epoch. The amplitude of the spectral features, both at low and high frequencies, directly depends 
on the value of the y-parameter, which describes the intrinsic CMB spectral distortion resulting from the energy release. Also the 
time-dependence of the injection process play an important role, for example leading to non-trivial shifts of the quasi-periodic pattern 
at low frequencies along the frequency axis. The existence of these narrow spectral features would open an unique way to separate 
y-distortions due to pre-recombinational (1400 < z < 5 x 10 4 ) energy release from those arising in the post-recombinational era at 
redshifts z < 800. 
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1. Introduction Also in the post-recombinational epoch (z < 

The measurements with the Cobe/Firas instrument have proven 
that the spectrum of the cosmic micro wave background (C MB) 
is very close to a perfect blackbody dFixsen et alJ |1996|) with 
thermodynamic temperature T q = 2.725 + 0.001 K( Mathe r et alJ 
■ Il999t iFixsen & Mather |2002|) . However, from the theoretical 
point of view deviations of the CMB spectrum from the one 
of a pure blackbody are not only possible but even inevitable 
if, for example, energy was released in the early Universe (e.g. 
due to viscous damping of acoustic waves, or annihilating or de- 
caying particles). For very early energy release (5 x 10 4 < z < 
2 x 10 6 ) the resulting spectral dist ortion can be characterized as 
a Bose-Einstein yu-type distort ion dSunvaev & Zeldovich| [T970b; 
Illari onov & Svunvaevll975a llbl). while for energy release at low 
redshift s (z < 5 X 10 4 ) the distorti on is close to a y-type dis- 
tortion dZeldovich & Sunvaevlll969l) . The current best observa- 
tional limits on th ese types of disto rtions are |y] < 1.5 x 10~ 5 and 
\p\ < 9.0 x 10" 5 dFixsen et al.1 1 19961) . Due to the rapid techno- 
logical development, improvements of these limits by a factor of 
~ 5 in principle could have been possible already several years 
ago dFixsen & Mathe rf2002l) . and recently some efforts are made 
to determine the absolute value of the CMB brightness tempera - 
ture at low frequencies using Arcade (Kog ut et alj|2004l l2006). 
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type spectral distortions due to different physical mechanisms 
should be produced. Performing measurements of the aver- 
age CMB spectrum (e.g. with wide-angle horns or like it was 
done with Cobe/Firas) all clusters of gala xies, hosting hot inter- 
galact ic gas, due to the thermal SZ-effect dSunvaev & Zeldovichl 
Il972bl) are contributing to the integral value of the observed 
y-parameter. S imilarly supernova remnants at high redshifts 
(lOh et al.ll2003l) . or shock waves arising due to l arge-scale struc- 
ture f ormation flSunvaev & Zeldovichl Il972at ICen & Os triker 
119991: iMiniati et al.1 12000|) should contribute to the overall y- 
parameter. For its value today we only have the upper limit by 
Cobe/Firas, and lower limits derived by estimating the total con- 
tribution of all clusters in the Universe dMarkevitch et"ai1 1 1 991; 
Ida Silva et al.ll2000t iRoncarelli et al.ll2.007h . These lower limits 
are exceeding y ~ 10~ 6 , and it is still possible that the contribu- 
tions to the total value of y due to early energy release are compa- 
rable or exceeding those coming from the low redshift Universe. 

Several detailed analytical and numerical studies for 
various energy injection histori es and mechanisms c an be 
found in the litera ture (e.g. Zeldovich & Sunyaev 1969; 
Sunyaev & Zeldovichl Il970bldlal ; llllarionov & Svunvaeyl 



Zeldovich et al.l 1 1 972t iChan et al.1 
Daly 1991; Burigana et al. 
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Danese & de Zottil 119821; iDalvl 119911: iBurigana et al.1 Il991bll; 
1995; Hu & Silk 1993a,b; Hu et al. 1994; Salvaterra & Burigan 
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Two very important conclusions can be drawn from these all 
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studies: (i) the arising spectral distortions are always very broad 
and practically featureless, and (ii) due to the absence of narrow 
spectral features, distinguishing different injection histories is 
extremely difficult. This implies that if one would find a y-type 
spectral distortion in the average CMB spectrum, then it is prac- 
tically impossible to say if the energy injection occurred just 
before, during or after the epoch of cosmological recombination. 
In this paper we show that the pre-recombinational emission 
within the bound-bound and free-bound transition of atomic 
hydrogen and helium should leave multiple narrow features 
(Av/v ~ 10 - 30%) in the CMB spectrum, that might become 
observable at cm, dm and sub-mm wavelength (see Sect. [5]). 
This could in principle open a way to directly distinguish pre- 
and post-recombinational y-distortions and even to shed light on 
the time-dependence of the energy injection process. 

At redshifts well before the epoch of He in — > He n- 
recombination (z > 8000) the total number of CMB pho- 
tons is not affected by atomic transitions if the intrinsic CMB 
spectrum is given by a pure blackbody. This is because the 
atomic emission and absorption processes in full thermody- 
namic equilibrium balance each other. However, a lower red- 
shifts (z < 8000), due to the expansion of the Universe, 
the medium became sufficiently cold to allow the formation of 
neutral atoms. The transition to the neutral state is associated 
with the release of several addition al photons per baryon (e. g. 
~ 5 photons per hydrogen atom dChluba & Sunvaevl 12006)), 
even within a pure blac kbody ambient CMB radiation field - 
Refining early estimates ([Zeldovich et al. 19681 : Pee bles 1968; 
lDubrovicHl975tlDubrovich & Stolvarovlll995Lll997h . the spec- 
tral distortions arising during hydrogen recombination (800 < 
Z < 1800), He II -> He i-recombination (1600 < z < 3000), and 
He in -> He n-recombination (4500 < z < 7000) within a pure 
blackbod y ambient radiation field have been recently computed 
in detail dRubino-Martm et all 120061: IChluba & Sunvaevl 120061 
IChluba et alJl2007aHRubino-M artin et alj|2007l) . It was also em- 
phasized that measuring these distortions in principle may open 
another independent way to determine the temperature of the 
CMB monopole, the specific entropy of the Universe, and the 
primordial helium abundance, well before the first appearance 
of stars (e.g.lSunvaev & Chluball2007l:IChluba & Sunvaevl2008t 
ISunvaev & Chluball2008l) ~ 

If on the other hand the intrinsic CMB spectrum devi- 
ates from a pure blackbody, full equilibrium is perturbed, 
and the small imbalance between emission and absorption in 
atomic transitions can lead to a net change of the number 
of photons, even prior to the epoch of recombination, in par- 
ticular owing to loops starti ng and ending in the continuum 
dLvubarskv & Sunvaevl [l983). These loops are trying to dimin- 
ish the maximal spectral distortions and are producing several 
new photons per absorbed one. In this paper we want to demon- 
strate how the cosmological recombination spectrum is affected 
if one allows for an intrinsic y-type CMB spectral distortion. 
We investigate the cases of single instantaneous energy injection 
(e.g. due to phase transitions) and for long-lived decaying par- 
ticles. There is no principle difficulty in performing the calcula- 
tions for more general injection histories, also including ^f-type 
distortions, if necessary. However, this still requires a slightly 
more detailed study, which will be left for a future paper. 

In Sect. [2] we provide a short overview regarding the ther- 
malization of CMB spectral distortions after early energy re- 
lease, and provide formulae which we then use in our compu- 
tations to describe y-type distortions. In Sect. [3] we give explicit 
expressions for the net bound-bound and free-bound rates in a 
distorted ambient radiation field. We then derive some estimates 



for the expected contributions to the pre-recombinational signals 
coming from primordial helium in Sect. [4] Our main results are 
presented in Sect. [3] where we first start by discussing a few 
simple cases (Sect. [57X1 and 15. 2t in order to gain some level of 
understanding. We support our numerical computations by sev- 
eral analytic considerations in Sect 15.1.11 and Appendix [B] In 
Sect. 15.31 we then discuss the results for our 25 shell computa- 
tions of hydrogen and He n. First we consider the dependence 
of the spectral distortions on the value of y (Sect. 15.3.11 ), where 
Fig.[5]and[6]play the main role. Then in Sects. |5372l and |5~3~3l we 
investigate the dependence of the spectral distortions on the in- 
jection redshift and history, where we are particularly interested 
in the low frequency variability of the signal (see Fig.l9land[TTli. 
We conclude in Sect. [6] 



2. CMB spectral distortions after energy release 

After any energy release in the Universe, the thermodynamic 
equilibrium between matter and radiation in general will be per- 
turbed, and in particular, the distribution of photons will de- 
viate from the one of a pure blackbody. The combi ned action 
of Compton scattering, double Compton emissiorQ (Lightman 
11981b lThornelll98ll: IChluba et alj|2007bl) . and Bremsstrahlung 
will attempt to restore full equilibrium, but, depending on the 
injection redshift, might not fully succ e ed. Usi ng t he approxi- 
mate formulae given in Buri eanaetall d!991bl) and[Hu & Silk 
( 1993a|), f or the parameters within the concordanc e cosmologi- 
cal model (Spergel et al. 2003; Be nnett et alj|2003l) . one can dis- 
tinguish between the following cases for the residual CMB spec- 
tral distortions arising from a single energy injection, 6p 7 /p y <sc 
1, at heating redshift Zh- 

(I) Zh < Z y ~ 6.3 x 10 3 : Compton scattering is not able 
to establish full kinetic equilibrium of the photon distribu- 
tion with the electrons. Photon producing processes (mainly 
Bremsstrahlung) can only restore a Planckian spectrum at 
very low frequencies. Heating result s in a Compton y- 
distortion (Zeld ovich & SunvaevH l969) at high frequencies, 
like in the case of the thermal SZ effect, with y-parameter 

(II) Zy < Zh < Zfi ~ 2.9 x 10 5 : Compton scattering can establish 
partial kinetic equilibrium of the photon distribution with 
the electrons. Photons that are produced at low frequencies 
(mainly due to Bremsstrahlung) diminish the spectral dis- 
tortion close to their initial frequency, but cannot upscatter 
strongly. The deviations from a blackbody represent a mix- 
ture between a y-distortion and a yu-distortion. 

(III) Zfi < Zh < z± ~ 2 x 10 6 : Compton scattering can estab- 
lish full kinetic equilibrium of the photon distribution with 
the electrons after a very short time. Low frequency pho- 
tons (mainly due to double Compton emission) upscatter 
and slowly reduce the spectral distortion at high frequen- 
cies The deviations from a blackbody can be described as a 
Bose-Einstein distribution with frequency-dependent chem- 
ical potential, which is constant at high and vanishes at low 
frequencies. 

(IV) Zth < Zh^ Both Compton scattering and photon production 
processes are extremely efficient and restore practically any 
spectral distortion arising after heating, eventually yielding 



1 Due to the huge excess in the number of photons over baryons 
(N r /Nb ~ 1.6 x 10 9 ), the double Compton process is the dominant 
source of new photons at redshifts zac <: 3 x 10 5 - 4 x 10 5 , while at 
z ^ Zdc Bremsstrahlung is more important. 
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a pure blackbody spectrum with slightly higher temperature 
T y than before the energy release. 

For the case (I) and (III) it is possible to approximate the dis- 
torted radiation spectrum analytically. There is no principle diffi- 
culty in numerically computing the time-dep endent solution for 
the radiation field after release of energy (e.g. lHu & SiM i993a) 
for more general cases, if necessary. However, here we are par- 
ticularly interested in demonstrating the main difference in the 
additional radiation appearing due to atomic transitions in hy- 
drogen and helium before the actual epoch of recombination, 
and for late energy release, in which a y-type spectral distortion 
is formed. We will therefore only distinguish between case (I) 
and (III), and assume that the transition between these two cases 
occurs at z^ v w z^l^^fl ~ 5.1 x 10 4 , i.e. the redshift at which 
the energy exchange time scale equals the expa nsion time scale 
dSunvaev & Zeldovich|[T980l:lHu & Silk|[T993ah . The effects due 
to intrinsic //-type spectral distortions will be left for a future 
work, so that below we will restrict ourselves to energy injection 
at redshift z < 50000. 



2. 1 . Compton y-distortion 

For energy release at low redshifts the Compton process is no 
longer able to establish kinetic equilibrium. If the temperature 
of the radiation is smaller than the temperature of the electrons 
photons are upscattered. For photons which are initially dis- 
tributed according to a blackbody spectrum with temperature T y , 
the efficiency of this process is determined by the Compton y- 
parameter, 



' k B (T, - T y ) 



y = J ^" c / : ' ,V,<r-,-d/. 



(1) 



where cr T is the Thomson cross section, dl = c dt, A e the elec- 
tron number density, and T e the electron temperature. For y «: 1, 
the resulting intrinsic distortion in the photo n occupation number 
of the CMB is approximately given by dZeldovich & Sunvaevl 
1969) 



Arty = y 



(e x - l) 2 



e x + 1 

x 4 

e x - 1 



(2) 



Here x = hv/k^Ty is the dimensionless frequency. 

For computational reasons it is convenient to introduce the 
frequency dependent chemical potential resulting from a y- 
distortion, which can be obtained with 



yu(x) = In 



1 + H, 



'y 



-y x 



e x + 1 
e x - 1 



(3) 



Here n p i(x) = l/[e x 



1] is the Planckian occupation number, 
and n 7 = « p i + An y . For x — > and }« 1 one finds p.(x) » 2xy, 
and for x » 1 one has p(x) as -ln[l + yx 2 ], or p(x) m -yx 2 
for 1 <K x «: -Jl/y. Comparing with a blackbody spectrum of 
temperature T y , for y > there is a deficit of photons at low 
frequencies, while there is an excess at high frequencies. In par- 
ticular, the spectral distortion changes sign at x y ~ 3.8. 



2.1.1. Compton y-distortion from decaying particles 

If all the energy is released at a single redshift, z, < z^ y ~ 50000, 
then after a very short time a y-type distortion is formed, where 
the y-parameter is approximately given by y ~ j 6p y /p y . 



However, when the energy release is due to decaying un- 
stable particles, which have sufficiently long life-times, f x , then 
the CMB spectral distortion will built up as a function of red- 
shift. In this case the fractional energy injection rate is given by 



Sf)ylPy 



a -m/tx 



/(l +z), so that the time-dependent y-parameter 



can be computed as 



y(z) = yo x 



>\2 



J z °° dz'e-Wx/Hiz'Xl+z') 
dz'e-Wx/H(z')(l +z') 2 ' 



(4) 



where yo - \ Sp y /p y is related to the total energy release, and 
H(z) is the Hubble expansion factor. Note that y(z) is a rather 
steep function of redshift, which strongly rises around the red- 
shift, zx, at which t(z) = fx- 



3. Atomic transitions in a distorted ambient CMB 
radiation field 

3. 1 . Bound-bound transitions 

Using the occupation number of photons, n y - l/[e x+fl — 1], with 
frequency-dependent chemical potential p{x), one can express 
the net rate connecting two bound atomic states i and j in the 
convenient form 



ARy = Pij 



AijNte-^i 

eXijWi - 1 



1 _ ii ^ e -\*fWi\ 
8j N, 



(5) 



where pij is the Sobolev-escape probability, Ay is the Einstein- 
A-coefficient of the transition i — > j, N, and g t are the popu- 
lation and statistical weight of the upper and Nj and gj of the 
lower hydrogen level, respectively. Furthermore we have intro- 
duced the dimensionless frequency Xu = hvij/kTo(l + z) of the 
transition, where 7p = 2 .725 K is the present CMB temperature 
(|Fixsen & Matherl l2002). and p. u = p.(x u ). 



3.2. Free-bound transitions 

For the free-bound transitions from the continuum to the bound 
atomic states i one has 



ARci^NeNcai-NiPi 



(6) 



where A^ c in the case of hydrogen is the number density of free 
protons, Af p , and the number density for Hem nuclei, A^Hem, in 
the case of helium. The recombination coefficient, a,, and pho- 
toionization coefficient, are given by the integrals 



a; = — fi(T e ) 



v 2 cr,-(v) e *+M*)+(%-*)/p 



8tt f°° vV,(v) 



dv (7a) 
(7b) 



Here x, c = hvi C /kT y is the dimensionless ionization frequency, 
p = TJTy is the ration of the photon and electron tempera- 
ture, cr, is the photoionization cross-section for the level i, and 

cm 3 . In full ther- 



•3/2. 



/to = f [sSht]' *|4.14xl0- 16 ^ 
modynamic equilibrium the photon distribution is given by a 
blackbody with T y = T e . As expected, in this case one finds 
from Eq. © that of = fi(T e ) e hv ^ kT 'pf. 
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4. Expected contributions from helium 



4.2. Contributions due to He i 



The number of helium nuclei is only ~ 8 % relative to the number 
of hydrogen atoms in the Universe. Compared to the radiation 
coming from hydrogen one therefore naively expects a small ad- 
dition of photons due to atomic transitions in helium. However, 
at given frequency the photons due to He n have been released 
at about Z 2 = 4 times higher redshifts than for hydrogen, so 
that both the number density of particles and temperature of the 
medium was higher. In addition the expansion of the Universe 
was faster. As we will show below, these circumstances make 
the contributions from helium comparable to those from hydro- 
gen, where He n plays a much more important role than He i. 

4.1. Contributions due to Hen 

The speed at which atomic loops can be passed through is deter- 
mined by the effective recombination rate to a given level i, since 
the bound-bound rates are always much faster. In order to esti- 
mate the contributions to the CMB spectral distortion by He n, 
we compute the change in the population of level i due to direct 
recombinations to that level over a very short time interval Af, 
i.e. ANi ~ N e N c af eu At. 

Because all the bound-bound transition rates in Hen are 16 
times larger than for hydrogen, the relative importance of the 
different channels to lower states should remain the same as in 
hydrogerfl Therefore one can assume that the relative number of 
photons, fij, emitted in the transition i — > j per additional elec- 
tron on the level i is like for hydrogen at 4 times lower redshifts. 

If we want to know how many of the emitted photons are 
observed in a fixed frequency interval Av today (z bs = 0) we 
also have to consider that at higher redshift the expansion of the 
Universe is faster. Hence the redshifting of photons through a 
given interval Av is accomplished in a shorter time interval. For 
a given transition, these are related by Af = „!t z Av. Then the 

change in the number of photons due to emission in the transition 
; — > j today should be proportional to 



a aw \ fijiZem) ANi(Zem) (1 + 

A/V y (V;;) Av, 

71 //(Zem)(l+Zem) 3 Vy 



(8) 



where z em is the redshift of emission, and the change of the vol- 
ume element due to the expansion of the Universe is taken into 
account by the factor of (1 + z em ) 3 - This now has to be compared 
with the corresponding change in the number of photons emitted 
in the same transition by hydrogen, but at 4 times lower emission 
redshift. 

For hydrogenic atoms with change Z the recombination rate, 
including stimulated r ecombination within the ambient CMB 
blackbody, scales like (Kaplan & Pike lnerll 19701) 



7-3/2 



f 

Jin 



hvi/kT 



dx 

72 



z 2 

7-1/2 



(9) 



where v, is the ionization frequency of the level i, T is the 
temperature of the plasma. It was assumed that /iv,- <K kT. 
Therefore one finds aP en (47')/arj II (7') ~ 2. Assuming radia- 
tion domination one also has H(z)/H(4z) ~ 1/16. Hence, we 
find AAr» eII (v^ eII ,4zem)/A/V^ iI (v|] I ,Zem) ~ 8%x4 3 2/16 ~ 64%. 
Note that (1 + ZemVvJ] 1 = (1 + 4z em )/v^ en . Prior to the epoch of 
He in — > He II recombination the release of photons by helium 
is amplified by a factor of ~ 8 ! 



In the case of neutral helium, the highly excited levels are ba- 
sically hydrogenic. Therefore one does not expect any ampli- 
fication of the emission within loops prior to its recombina- 
tion epoch. Furthermore, the total period during which neu- 
tral helium can contribute significantly is limited to the redshift 
range starting at the end of He in — > He n recombination, say 
1600 < z < 6000. Therefore, neutral helium typically will not 
be active over a large range of redshifts. 

Still there could be some interesting features appearing in 
connection with the fine-structure transitions, which even within 
the standard computations lea d to strong negative features in the 
He i recombination spectrum (Rubino-Martin et al. 2007). Also 
the spectrum of neutral helium, especially a high frequencies, is 
more complicated than for hydrogenic atoms, so that some non- 
trivial features might arise. We leave this problem for some fu- 
ture work, and focus on the contributions of hydrogen and He n. 

5. Results for intrinsic v-type CMB distortions 

Here we discuss the results for the changes in the recombina- 
tion spectra of hydrogen and He n for different values of the y- 
parameter. Some of the computational details and the formula- 
tion of the problem can be found in Appendix lAl 

5.1. The 2 shell atom 

In order to understand the properties of the solution and also to 
check the correctness of our computations we first considered 
the results including only a small number of shells. If we take 
2 shells into account, we are only dealing with a few atomic 
transitions, namely the Lyman- and B aimer-continuum, and the 
Lyman-ff line. In addition, one expects that during the recom- 
bination epoch of the considered atomic species (here Hi or 
Hen) also the 2s-ls-two-photon decay channel will contribute, 
but very little before that time. 

In Fig. [TJwe show the spectral distortion, AI V , including 2 
shells into our computations for different transitions as a func- 
tion of redshiffl It was assumed that energy was released in a 
single injection at z; = 50000, leading to y = 10~ 5 . All shown 
curves were computed using the ^-function approxim ation for 
the intensity (for details see Rubin o-Martfn et al.l2006l) . This ap- 
proximation is not sufficient when one is interested in computing 
the spectral distortions in frequency space. 

Prior to the recombination epoch of the considered species 
one can find pre-recombinational emission and absorption in the 
Lyman- and B aimer-continuum, and the Lyman-c line, which 
would be completely absent for y = 0. As expected, during in 
the pre-recombinational epochs the 2s-ls-two-photon transition 
is not important. This is because the 2s- Is transition is simply 
unable to compete with the ~ 10 8 times faster Lyman-a transi- 
tion while it is still optically thin. 

Summing the spectral distortions due to the continuua, one 
finds cancellation of the redshift-dependent emission at a level 
close to our numerical accuracy (relative accuracy < 10~ 4 for 
the spectrum). This is expected because of electron number con- 
servation: in the pre-recombinational epoch the overall ioniza- 
tion state of the plasma is not affected significantly by the small 



2 Even the factors due to stimulated emission in the ambient black- 
body radiation field are the same! 



3 This is a convenient representation of the spectrum, when one is 
interested in the time-dependence of the photon release, rather than the 
observed spectral distortion in frequency space. To obtain the later, one 
simply has to plot the presented curves as a function of v = Vjj/(1 + z), 
where is the restframe frequency of the considered transition. 
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Fig. 1. Spectral distortion, A/ v , including 2 shells into our com- 
putations for different transitions as a function of redshift and for 
y - 10~ 5 . All shown curves were computed using the ^-function 
approximation for the intensity. The upper panel shows the re- 
sults for hydrogen, the lower those for He n. In both cases the 
analytic approximation for the Lyman-*? line based on Eq. < |B.2K 
and including the escape probabilities in the Lyman-cr line and 
Lyman-continuum, are also shown. 
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Fig. 2. Analytic representation of the pre-recombinational Hi 
Lyman-o' spectral distortions for the 2 shell atom and y - 10 -5 . 
See text for explanations. 



deviations of the background radiation from full equilibrium. 
Therefore all electrons that enter an atomic species will leave 
it again, in general via another route to the continuum. This im- 
plies that 2,- AR C( - = 0, which is a general property of the solution 
in the pre-recombinational epoch. 

If we look at the Lyman- and Balmer-continuum in the case 
of hydrogen we can see that at redshifts z > 30000 electrons are 
entering via the Lyman-continuum, and leaving via the Balmer- 
continuum, while in the redshift range 2000 < z < 30000 the 
opposite is true. As expected, for z > 2000 the Lyman-a transi- 
tion closely follows the Balmer-continuum, since every electron 
that enters the 2p-state and then reaches the ground level, had 
to pass through the Lyman-a transition. Using the analytic so- 
lution for the Lyman-c line as given in the Appendix iBl we find 
excellent agreement with the numerical results until the real re- 
combination epoch is entered at z < 2000. 

In the case of He n for the considered range of redshifts the 
pre-recombinational emission (z > 7000) is always generated in 
the loop c — > 2p — > Is — > c. Again we find excellent agree- 
ment with the analytic solution for the Lyman-o. line. Note that 
for He ii the total emission in the pre-recombinational epoch is 
much larger than in the recombination epoch at z ~ 6000 (see 
discussion in Sect. The height of the maximum is even com- 
parable with the Hi Lyman-a' line. 

As one can see from Fig. Q] at high redshifts all transitions 
become weaker. This is due to the fact that the restframe frequen- 
cies of all lines are in the Rayleigh- Jeans part of the CMB spec- 
trum, where the effective chemical potential of the y-distortion 
(see Sect. 12.11 ) is dropping lik^3ju(x) « 2xy. This implies that 
at higher redshift all transitions are more and more within a 
pure blackbody ambient radiation field. On the other hand the 
effective chemical potential increases towards lower redshift, so 
that also the strength of the transitions increases. However, at 
z < 3000 in the case of hydrogen, and z < 1 1000 for He n, the 
escape probability in the Lyman-continuum (see Appendix lA.il 
and Eq. ( |B.3t for quantitative estimates) starts to decrease sig- 
nificantly, so that the pre-recombinational transitions cease. The 
maximum in the pre-recombinational Lyman-ff line is formed 
due to this rather sharp transitions to the optically thick region 
in the Lyman-continuum (see also Sect. l5.1.Tl for more details). 

5.1 .1 . Analytic description of the pre-recombinational 
Lyman-Q- line 

One can understand the behavior of the solution for the spectral 
distortions in more detail using our analytic description of the 
Lyman-c line as given in AppendixlBl 

In Fig. [2] we show the comparison of different approxima- 
tions with the full numerical result. If we use the analytic ap- 
proximation based on Eq. SB. 2} , but do not include the escape 
probabilities in the H i Lyman-a' line and H i Lyman-continuum, 
then we obtain the dotted curve (quoted 'analytic la' in the fig- 
ure). The curve quoted 'analytic lb' also includes the escape 
probabilities as described in Appendix IB. 1.1 1 Comparing these 
curves shows that for the shape of the distortion at z < 3000 the 
escape probabilities are very important. However, although at 
this redshift the Sobolev optical depth in the H i Lyman-o' line is 
roughly 14 times larger than the optical depth in the H i Lyman- 
continuum, the derivation of Eq. dB- 8t shows that the Hi Lyman- 
a escape probability only plays a secondary role. 



4 Or more correctly fi(x) « 7.4 xy if one also takes into account the 
difference in the photon and electron temperature T e ss T y [l + 5Ay\ 



(see Appendix I A. 1 
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Fig. 3. Spectral distortion, A7 V , including 3 shells into our com- 
putations for different transitions as a function of redshift and 
for y = 10 . For all shown curves we used the (5-function ap- 
proximation to compute the intensity. The upper panel shows the 
results for hydrogen, the lower those for He rx. 



With the formulae in Appendix lB.1.21 the spectral distortion 
can be written in the form A/ V (z) = F(z) x A. If we use F ac- 
cording to Eq. ( IB. 10| > and A ss ^21 + A%>c - A'isc, as derived in 
Eq. ( lB.14bl >. then we obtain the approximation quoted 'analytic 
IF. One can clearly see that this approximation represents the 
global behavior, but is fails to explain the Lyman-or absorption 
at z > 30000. In fact within this approximation the Lyman-o' line 
should always be in emission, even at very high redshifts, since 
there A * 3yxi sc /32 > 0. 

If we also take into account higher order terms for the line 
imbalance A according to Eq. ( IB. 16b then we obtain the curve 
quoted 'analytic III', which is already very close to the full solu- 
tion and also reproduces the high redshift behavior, but starting 
at slightly higher redshift (z ~ 40000 instead of z ~ 3 0000) 
This is largely due to the approximations of the integrals ( lB.15t 
over the photoionization cross-sections (in particular M_i). Still 
if one evaluates these integrals more accurately, one does not re- 
cover the full solution, since the free-bound Gaunt-factors were 
neglected. 

5.2. The 3-shell atom 

If one takes 3 shells into account, the situation becomes a bit 
more complicated, since more loops connecting to the contin- 



Fig. 4. Sketch of the main atomic loops for hydrogen and He 11 
when including 3 shells. The left panel shows the loops for 
transitions that are terminating in the Lyman-continuum. The 
right panel shows the case, when the Lyman-continuum is com- 
pletely blocked, and unbalanced transitions are terminating in 
the B aimer-continuum instead. 



uum are possible. Looking at Fig.[3]again we find that the sum 
over all transition in the continuua vanishes at redshifts prior 
to the actual recombination epoch of the considered species. At 
z < 3000 in the case of hydrogen, and z < 11000 for Hen, 
the escape probability in the Lyman-continuum becomes small. 
For 2 shells this fact stopped the pre-recombinational emission 
until the actual recombination epoch of the considered atomic 
species was entered (see Fig. [T}. However, for 3 shells electrons 
can now start to leave the ls-level via the Lyman-yS transition, 
and then reach the continuum through the Balmer-continuum. 
For both hydrogen and He n one can also see that the emission in 
the Lyman-a line stops completely, once the Lyman-continuum 
is fully blocked. In this situation only the loop c — > 3 — > 2 — > c 
via the Balmer-continuum is working. Only when the main re- 
combinational epoch of the considered species is entered, the 
Lyman-or line is reactivated. 

In Fig. |4]we sketch the main atomic loops in hydrogen and 
Hen when including 3 shells. For y = 10 , in the case of hy- 
drogen the illustrated Lyman-continuum loops work in the red- 
shift range 2000 < z < 30000, while the Balmer-continuum 
loop works for 1600 < z < 2000. In the case of Hen one finds 
8000 < z < 1.2 x 10 5 and 6200 < z < 8000 for the Lyman- 
and Balmer-continuum loops, respectively. It is clear that in ev- 
ery closed loop one energetic photons is destroyed and at least 
two photons are generated at lower frequencies. Including more 
shells will open the possibility to generate more photons per 
loop, simply because electrons can enter through highly excited 
levels and then preferentially cascade down to the lowest shells 
via several intermediate levels, leaving the atomic species taking 
the fastest available route back to the continuum. Below we will 
discuss this situation in more detail (see Sect. |5.3.TT >. 

Figures[T]and|3]both show that the pre-recombinational lines 
are emitted in a typical redshift range Az/z ~ 1, while the signals 
from the considered recombinational epoch are released within 
Az/z ~ 0.1 - 0.2. For the pre-recombinational signal the ex- 
pected line-width is Av/v ~ 0.6 - 0.7. However, the overlap of 
several lines, especially at frequencies where emission and ab- 
sorption features nearly coincide, and the asymmetry of the pre- 
recombinational line profiles, still leads to more narrow spectral 
features with Av/v ~ 0.1 - 0.3 (see Sect.ES Fig.©. 

It is also important to mention that in all cases the actual 
recombination epoch is not affected significantly by the small y- 
distortion in the ambient photon field. There the deviations from 
Saha-equilibrium because of the recombination dynamics dom- 
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Fig. 5. Contributions to the H i (left panels) and He n (right panels) recombination spectrum for different values of the initial y- 
parameter. Energy injection was assumed to occur at z, ■ - 4 x 10 4 . In each column the upper panel shows the bound-bound signal, 
the middle the free-bound signal, and the lower panel the sum of both. The thin red lines represent the overall negative parts of the 
signals. 



inate over those directly related to the spectral distortion, and in 
particular the changes in the ionization history are tiny. 



5.3. The 25-shell atom 

In this Section we discuss the results for our 25-shell compu- 
tations. Given the large amount of transitions, it is better to di- 
rectly look at the spectral distortion as a function of frequency. 
However, following the approach of Sect 15. II we have checked 
that the basic properties of the first few lines and continuua as a 
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function of redshift do not change qualitatively in comparison to 
the previous cases. In particular, only the Lyman- and Balmer- 
continuum become strongly negative, again in similar redshifts 
ranges as for 2 or 3 shells. In absorption the other continuua play 
no important role. 

Note that for the spectral distortion now the impact of elec- 
tron scattering has to be considered and the free-bound contribu- 



tion must be computed using the full differential photoionization 
cross section (see Appendix |A.4| ). 

5.3.1 . Dependence of the distortion on the value of y 

In Fig. [5] we show the contributions to the recombination spec- 
trum for different values of the initial y-parameter. In addition, 
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Fig. [6] shows some of the main contributions to the total hydro- 
gen and He n spectral distortion in more detail. 



Bound-bound transitions 

Focusing first on the contributions due to bound-bound transi- 
tions, one can see that the standard recombination signal due to 
hydrogen is not very strongly affected when y = 10~ 7 , whereas 
the helium signal already changes notably. Increasing the value 
of y in both cases leads to an increase in the overall amplitude of 
the distortion at low frequencies, and a large rise in the emission 
and variability at v > 100 GHz. For Hi at low frequencies the 
level of the signal changes roughly 5 times when increasing the 
value of y from to 10~ 5 , while for Hen the increase is even 
about a factor of 40. This shows that in the pre-recombinational 
epoch He n indeed behaves similar to hydrogen, but with an am- 
plification ~ 8 (see Sect. |4j>. 

At high frequencies, a strong emission-absorption feature 
appears in the range v ~ 500 GHz - 1600 GHz, which is com- 
pletely absent for y — 0. With y = 10~ 5 from peak to peak 
this feature exceeds the normal Lyman-o' distortion (close to 
v ~ 1750GHz for Hi and v ~ 1680GHz for Hen) by a fac- 
tor of ~ 5 for H i, and about 30 times for He n. The absorption 
part is mainly due to the pre-recombinational Lyman-/?, -y, and - 
6 transition, while the emission part is dominated by the pre- 
recombinational emission in the Lyman-** line (see also Fig. [6] 
for more detail). 

It is important to note that in the case of He n most of the re- 
combinational Lyman-o' emission (v ~ 1750 GHz) is completely 
wiped out by the pre-recombinational absorption in the higher 
Lyman-series, while for H i only a very small part of the Lyman- 
a low frequency wing is affected. This is possible only because 
the pre-recombinational emission in the He n Lyman-series is so 
strongly enhanced, as compared to the signal produced during 
the recombinational epoch. 



Free-bound transitions 

Now looking at the free-bound contributions, one can again see 
that the hydrogen signal changes much less with increasing value 
of y than in the case of He n. In both cases the variability of the 
free-bound signal decreases at low frequencies, while at high fre- 
quencies a strong and broad absorption feature appears, which 
is mainly due to the Lyman-continuum. For y = 10~ 5 this ab- 
sorption feature even completely erases the Balmer-continuum 
contribution appearing during the actual recombination epoch of 
the considered species. It is 2 times stronger than the H i Lyman- 
a line from the recombination epoch, and in the case of He n it 
exceeds the normal He n Lyman-*? line by more than one order 
of magnitude. 

However, except for the absorption feature at high frequen- 
cies the free-bound contribution becomes practically feature- 
less when reaching y = 10 . This is due to the strong over- 
lap of different lines from the high redshift part, since the char- 
acteristic width of the recombinational emission increases like 
Av/v ~ kT e /hvj C (see middle panels in Fig.|5j. In addition, the 
photons are released in a more broad range of redshifts (see 
Sect. 15. II and 15. 21 . also leading to a lowering of the contrast of 
the spectral features from the recombinational epoch. 



Table 1. Approximate number of photons and loops per nucleus 
for Zi ■ - 4 x 10 4 and different values of y. 



y = y = IP' 7 y = IP' 6 y = 1(T T 

AF^ Ox N ti -0.028 x N H -0.28 x N H -2.74 x N H 

N$_ a 0.42 xN» 0.47 xA^h 0.92 x N H 5.37 x N H 

Aff b ' 2.49 xN H 2.59 x N H 3.45 x N H 12.04 x N tt 

A^ o ' p 0xN H 0.032 xA'h 0.32 xN H 3.24 x~/Vh~ 

OxA^He -0.27 x A'He -2.67 x N Hc -26.6 x N Hc 

A7gL° 0.55 x ^ 1.06 x A'He 5.68 x N He 51.8 x^ 

Affi" 2.48 x A'He 3.49 x A'He 12.5 x N Hs 103 x N Hs 

WE 1 0^N^ 0.28 x A'He 2.91 x A^e 30.1 xA^e 



Total distortion 

Also in the total spectra (see lower panels in Fig. |5]l one can 
clearly see a strong absorption feature at high frequencies, which 
is mainly associated with the Lyman-continuum and Lyman- 
series for n > 2 (see Fig. [6] also). For y = 10 5 , in the case 
of hydrogen it exceeds the Lyman-a line from recombination 
(v ~ 1750 GHz) by a factor of ~ 4 at v ~ 1250 GHz, while for 
He ii it is even ~ 20 times stronger, reaching ~ 80% of the cor- 
responding hydrogen feature. Checking the level of emission at 
low frequencies, as expected (see Sect.[4j>, one can find that He n 
indeed contributes about 2/3 to the total level of emission. 

As illustrated in the upper panels of Fig. |6l the emission- 
absorption feature at high frequencies is due to the overlap of the 
pre-recombinational Lyman-a line (emission), and the combina- 
tion of the higher pre-recombinational Lyman-series and Lyman- 
continuum (absorption). At intermediate frequencies (middle 
panels), the main spectral features are due to the Balmer-a, pre- 
recombinational Balmer-series from n > 3 and the Paschen-ar 
transition, with some additional broad contributions to the over- 
all amplitude of the bump coming from higher continuua. 

The lower panels of Fig. [6] show, the separate contributions 
to the bound-bound series for the 10th shell as an example. One 
can notice that in the case of hydrogen the recombinational and 
pre-recombinational emission have similar amplitude, while for 
He ii the pre-recombinational signal is more than one order of 
magnitude larger (see Sect. |4]for explanation). In both cases the 
pre-recombinational emission is much broader than the recom- 
binational signal, again mainly due to the time-dependence of 
the photon emission process (see Sect. 15. 1 | and |5.2t , but to some 
extent also because of electron scattering. 

Number of photons and loops 

Using the free-bound spectrum, one can also compute the total 
number of loops, A/ioop> that were involved into the production 
of photons. This is possible, since only during the recombination 
epoch electrons will terminate in the ls-state. Therefore the total 
number of photons emitted in the free-bound spectrum is very 
close to ~ \y per nucleus, while when looking at the positive or 
negative part of this contribution one should find ~ (A/i 00 p + 1)7 

and A/ioop7 per nucleus, respectively. 

In Table [TJ we give a few examples, also comparing with the 
number of photons emitted for y = 0. One can see that the num- 
ber of loops per nucleus scales roughly proportional to the values 
of y, i.e. ACp ~ 3-2 x [y/lQ- 5 ] and N^ 1 ~ 30 x [y/lQ- 5 ]. If 
one would consider a lower injection redshift the proportional- 
ity constant should decrease. Also when including more shells 
Wi 00 p should become larger. Furthermore, the number of loops 
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Fig. 7. H i (upper panel) and He n (lower panel) recombination 
spectra for different energy injection redshifts. 



per nucleus is about one order of magnitude larger for Hen 
than for hydrogen. As explained in Sect. |H this is due to the 
amplification of transitions in the case of hydrogenic helium at 
high redshifts. Comparing the number of photons absorbed in 
the Lyman-continuum, one can see that in practically all cases 
~ 90% of all loops are ending there. 

If we take the total number of photons per nucleus emitted in 
the bound-bound transitions and subtract the number of photons 
emitted for y — 0, we can estimate the loop-efficiency, ei 00p , or 
number of bound-bound photons generated per loop prior to the 
recombination epoch. For hydrogen one finds ei oop ~ 2.9 - 3.1, 
while for Hen one has ei oop ~ 3.3 - 3.6. Similarly one obtains 
an loop efficiency of ei oop ~ 1.7-1.8 for both the Hi and 
Hen Lyman-a lines. As expected these numbers are rather in- 
dependent of the value of y, since they should reflect an atomic 
property. They should also be rather independent of the injec- 
tion redshift, which mainly affects the total number of loops and 
thereby the total number of emitted photons. However, the loop 
efficiency should still increase when including more shells in the 
computation. 

5.3.2. Dependence of the distortion on the redshift of energy 
injection 

To understand how the pre-recombinational emission depends 
on the redshift at which the energy was released, in Fig. [7] we 




Fig. 8. Total Hi + He n recombination spectra for different en- 
ergy injection redshifts. The upper panel shows details of the 
spectrum at low, the lower at high frequencies. 



show a compilation of different cases for the total H i and He n 
signal. In Fig. [8] we also present the sum of both in more detail. 



Features at high frequencies 

For all shown cases the absolute changes in the curves are 
strongest at high frequencies (v > 100 GHz). Generally speak- 
ing, one can again find a rather broad bump at 100 GHz < 
v < 400 GHz, followed by an emission-absorption feature in 
the frequency range 500 GHz < v < 1600 GHz. In particular 
the strength and position of this emission-absorption feature de- 
pends strongly on the redshift of energy injection. 

Regarding the broad high-frequency signature, it is more im- 
portant that the variability is changing, rather than the increase 
in the overall amplitude. For example, in the frequency range 
100 GHz < v < 400 GHz the normal recombinational signal has 
~ 2 spectral features, while for injection at z, ~ 4000 roughly 
4 features are visible, which in this case only come from hydro- 
gen, since at z ~ 4000 He n is already completely recombined. 
Note that also neutral helium should add some signal, which was 
not included here. Nevertheless, we expect that this contribution 
is not strongly amplified like in the case of He n (see Sect. |4j, 
and hence should not add more than 10 - 20% to the total signal. 
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Fig. 9. Comparison of the variable component in the Hi + He n bound-bound and free-bound recombination spectra for single 
energy injection at different redshifts. In all cases the computations were performed including 25 shells and y = 10 5 . The blue 
dashed curve in all panel shows the variability in the normal H i + He n recombination spectrum (equivalent to energy injection 
below z ~ 800 or y = 0) for comparison. 



Variability at low frequencies 

Focusing on the spectral distortions at low frequencies, the over- 
all level of the distortion in general increases for higher red- 
shifts of energy release. However, there is also some change in 
the variability of the spectral distortions. In order to study this 
variability in more detail, in each case we performed a smooth 
spline fit of the total Hi + Hen recombination spectrum and 
then subtracted this curve from the total spectrum. The remain- 
ing modulation of the CMB intensity was then converted into 
variations of the CMB brightness temperature using the relation 
AT/To = hi I B v , where B v is the blackbody intensity with tem- 
perature 7q = 2.725 K. 

Figure|9]shows the results of this procedure for several cases. 
It is most striking that the amplitude of the variable component 
decreases with increasing energy injection redshift. This can be 
understood as follow: we have seen in Sect. 15.11 and 15.21 that 
for very early energy injection most of the pre-recombinational 
emission is expected to arise at z ~ 3000 for hydrogen, and 
z ~ 11000 for helium, (i.e. the redshifts at which the Lyman- 
continuum of the considered atomic species becomes optically 
thick) with a typical line-width Av/v ~ 1. In this case the total 



variability of the signal is mainly due to the non-trivial superpo- 
sition of many broad neighboring spectral features. Most impor- 
tantly, very little variability will be added by the high redshift 
wing of the pre-recombinational lines and in particular the be- 
ginning of the injection process. This is because (i) at high z the 
emission is much smaller (cf. Fig.Q]and[3]), and (ii) electron scat- 
tering broadens lines significantly, smoothing any broad feature 
even more (see Appendix |A.4| ). 

On the other hand, when the energy injection occurs at lower 
redshift this increases the variability of the signal because (i) 
electron scattering in the case of single momentary energy re- 
lease does not smooth the step-like feature due to the beginning 
of the injection process as strongly, and (ii) the total emission 
amplitude and hence the step-like feature increases (see Fig Q] 
and [3]). Once the injection occurs very close to the redshift at 
which the Lyman-continuum is becoming optically thick (see 
Sect. 15.11) . i.e. where the pre-recombinational emission has an 
extremum, this therefore should lead to a strong increase of the 
variability. On the other hand, for energy injection well before 
this epoch the atomic transitions lead to an increase in the over- 
all amplitude of the distortions rather than the variability. 
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Fig. 10. Comparison of the variable component in the standard 
(y — 0) H i + He n bound-bound and free-bound recombination 
spectrum for n max = 100 and 25. 



Indeed this can be also seen in Fig. [9] where for z, = 4000 
the variable component is ~ 3-8 times larger than the nor- 
mal recombinational signal, with a peak to peak amplitude of 
~ 50 - 70 nK instead of ~ 10 - 15 nK at frequencies around 
~ 1.5 GHz. Even for z, > 15000 the amplitude of the variable 
component is still 1.5-2 times larger than in the case of standard 
recombination, but it practically does not change anymore when 
going to higher injection redshifts. For z, ~ 11000 one expects 
a similarly strong increase in the variability as for Zi ~ 4000, 
but this time due to Hen. In addition to the change in ampli- 
tude of the variable component, in all cases the signal is shifted 
with respect to the normal recombinational signal. These shifts 
should also make it easier to distinguish the signatures from pre- 
recombinational energy release from those arising because of 
normal recombination. 

It is important to mention that the total amplitude of the 
variable component, should still increase when including mor e 
shells into the computation. As shown in lChluba et alj (12007a), 
for y = in particular the overall level of recombinational emis- 
sion at low frequencies strongly depends on the completeness of 
the atomic model. Similarly, the variability of the recombination 
spectrum changes. We illustrate this fact in Fig. [TUl where we 
compare the variability in the Hi + Hen recombination spec- 
trum for 25 sh ells (y = 0), with the one obtained in our 100 shel l 
computations dChluba et alj|2007al: iRubino-Martin et al.ll2007l) . 
As one can see, at low frequencies (v ~ 1-3 GHz) the amplitude 
of the variable component increases by more than a factor of 2 
when including 100 shells, reaching a peak to peak amplitude of 
~ 30 nK. This is due to the fact that for a more complete atomic 
model additional electrons are able to pass through a particular 
transition between highly excited states. 

5.3.3. Dependence of the distortion on the energy injection 
history 

Until now we have only considered cases for single momentary 
energy injection. However, physically this may not be very re- 
alistic, since most of the possible injection mechanisms release 
energy over a more broad range of redshifts. Also the discussion 
in Sect. 15.3.21 has shown that for single injection a large part of 
the variability can be attributed to the inset of the energy release. 
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Fig. 11. Comparison of the variable component in the H i + He n 
bound-bound and free-bound recombination spectra for single 
energy injection (black solid curves) and energy injection due 
to long-lived decaying particles with different lifetimes (red 
dashed-dotted curves). In all cases the computations were per- 
formed including 25 shells and a maximal y-parameter y - 10 . 



Therefore it is important to investigate the potential signatures 
of other injection mechanisms. 

For the signals under discussion long lived decaying particles 
are most interesting. In Sect. l2.1.T1 we have given some simpli- 
fied analytic description of this problem. In Fig.[TT]we show the 
variable component for the CMB spectral distortion due to the 
presence of hydrogen and He n at low frequencies, in the case 
of single injection and for energy release due to long-lived de- 
caying particles. It is clear that for zx = 4000 the variability 
is significantly smaller than in the case of single energy injec- 
tion. This is due to the fact that the inset of the atomic transi- 
tions is much more gradual than in the case of single injection. 
However one should mention that for energy injection due to 
decaying particles the effective y-parameter at z = 4000 is still 
only ~ 1 /3 of the maximal value, so that the level of variabil- 
ity cannot be directly compared with the case of single energy 
injection. Nevertheless, the structure of the variable component 
still depends non-trivially on the effective decay redshift, so that 
one in principle should be able to distinguish different injection 
scenarios. 
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Similarly, one could consider the case of annihilating parti- 
cles. However, here energy is effectively released at higher red- 
shiftsS and also within a much broader redshift interval. In this 
case, one has to follow the evolution of the CMB spectrum due 
to this heating mechanism from an initial /U-type distortion to a 
partial y-type distortion in more detail. Also one can expect that 
the redistribution of photons via electron scattering will become 
much more important (see Appendix lA.4l i. and that the free-free 
process will strongly alter the number of photons emitted via 
atomic transitions (see Appendix |A.51 >. In addition, collisional 
processes may become significant, in particular those leading to 
transitions among different bound-bound level, or to the contin- 
uum, since they are not associated with the emission of photons. 
This problem will be left for some future work. 

6. Discussion and conclusions 

In the previous Sections we have shown in detail how intrinsic y- 
type CMB spectral distortions modify the radiation released due 
to atomic transitions in primordial hydrogen and Hen at high 
redshifts. We presented the results of numerical computations 
including 25 atomic shells for both H i and He n, and discussed 
the contributions of several individual transitions in detail (e.g. 
see Fig. [6), also taking the broadening of lines due to electron 
scattering into account. As examples, we investigate the case of 
instantaneous energy release (Sect. |5.3.2t and exponential en- 
ergy release (Sect. |5.3.3t due to long lived decaying particles, 
separately. 

Our computations show that several additional photons are 
released during the pre-recombinational epoch, which in terms 
of number can strongly exceed those from the recombinational 
epoch (see Sect. 15.3. It . The number of loops per nucleus scales 
roughly proportional to the values of y, i.e. ~ 3.2x[y/10 -5 ] 

and Af^jp 1 ~ 30 x [y/10 ] for hydrogen and Hen respectively, 
where effectively about 3 photons per loop are emitted in the 
bound-bound transition. 

Due to the non-trivial overlap of broad neighboring pre- 
recombinational lines (from bound-bound and free-bound tran- 
sitions), rather narrow (Av/v ~ 0.1 -0.3) spectral features on top 
of a broad continuum appear, which both in shape and amplitude 
depend on the time-dependence of the energy injection process 
and the value of the intrinsic y-type CMB distortion. At high 
frequencies (v ~ 500 GHz - 1600 GHz) an emission-absorption 
feature is formed, which is completely absent for y = 0, and is 
mainly due to the superposition of pre-recombinational emission 
in the Lyman-or line, and the higher Lyman-series and Lyman 
continuum. Looking at Fig.Q~2]it becomes clear that this absorp- 
tion feature (close to v ~ 1400 GHz) in all shown cases even 
exceeds the intrinsic y-distortion. For y = 10 5 it even reaches 
~ 10% of the CMB blackbody intensity. Unfortunately, it ap- 
pears in the far Wien-tail of the CMB spectrum, where espe- 
cially the cosmic infra red background du e to dusty star- formin g 
galaxies is dominant dFixsen et al.lll998t lLagache et al.l 12005). 
Still one may hope to be able to extract such spectral feature in 
the future. 

One should stress that still all the discussed additional pre- 
recombinational spectral distortions are in general small in com- 
parison with the intrinsic y-distortion. As Fig. [T2l shows, the am- 
plitude of the additional distortions is typically well below 1% of 
the CMB y-distortion. However, at low frequencies (v ~ 1 GHz) 

5 This conclusion depends also on the temperature/energy depen- 
dence o f the annihilation cross s ection. We assumed i-wave annihilation 
(e.g. see lMcDonald et alj|2001l) . 




V [ GHz ] 

Fig. 12. Spectral distortions relative to the CMB blackbody spec- 
trum, B v . The thin blue curves show the absolute value of a 
y-type distortion. At low frequencies we indicate the expected 
level of emission when including more shells in our computa- 
tions. 



the additional distortions should reach ~ 10% of the intrinsic 
y-distortion, and a lower frequency may even exceed it. But in 
this context it is more important that due to the processes dis- 
cussed here new narrow spectral features appear, which have 
a unique variability (e.g. see Fig. that is even stronger than 
in the case of the recombinational lines from standard recombi- 
nation (y = 0). Such variability is very hard to mimic by any 
astrophysical foreground or systematic problem with the instru- 
ments. As em phasized earlier for the co smological recombina- 
tion spectrum ( Sun yaev & Ch luba 2007|, 120081) . this may allow 
us to measure them differentially, also making use of the fact that 
practically the same signal is coming from all directions on the 
sky. For intrinsic ^-distortions direct differential measurements 
are much harder, since it has so broad spectrum. Furthermore, 
as pointed out in the introduction, by measuring the narrow 
spectral features under discussion one could in principle distin- 
guish pre- and post-recombinational energy release, an observa- 
tion that cannot be achieved by directly measuring the average 
>'-distortion of the CMB. 

Above it was pointed out, that there is no principle difficulty 
in computing the spectral distortions due to pre-recombinational 
atomic transitions in H i and He n for more general energy in- 
jection histories, if necessary. In particular very early injection, 
involving //-type distortions, may be interesting, since stimulated 
emission could strongly amplify the emission a low frequen- 
cies, and hence the total number of emitted photons per atom. 
However, to treat this problem one has to follow the detailed evo- 
lution o f the CMB spectr al distortion due to the injection process 
(see e.g. lHu &" Silk 1993a). Also the effect of electron scattering 
on the distortions due to the pre-recombinational atomic transi- 
tions, in particular because of the recoil effect, has to be treated 
more rigorously. Simple estimates show (see Appendix |A.5t that 
at low frequencies also the modifications due to the free-free ab- 
sorption will become significant. Furthermore, one has to ac- 
count for collisional processes, since they should become im- 
portant at very high redshifts, even for shells with rather low n. 

An additional difficulty arises due to the fact that both at very 
low and at very high frequencies the back-reaction of the pre- 
recombinational distortion on the ambient radiation field may 
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not be negligible (see Fig. [T5V This could also still affect the 
details of the results presented here, but the main conclusions 
should not change. We defer all these problems to another paper. 
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Appendix A: Computational details 

A.1. Formulation of the problem 

Details on the formul ation of the problem for a p ure blackbody CMB radia- 
tion field are given in Rubino-Martfn et al. ( 2006) and references therein. We 
shall use the same notation here. The main difference with respect to the stan- 
dard recombination computation is due to the possibility of a non-blackbody 
ambient radiation field, which affects the net bound-bound and free-bound rates 
as explained in Sect. [5] Also the temperature of the electrons in general is no 
longer equal to the effective temperature of the photons, as we discuss below 
(see Sect. IA.31 . Since we consider only small intrinsic spectral distortions all 
the modification to the solution for the level populations are rather small, and 
most of the differences will appear only as pre-recombinational emission due to 
atomic transitions, but with practically no effect on the ionization history. 

One additional modificati on is r e lated to the Ly man-continuum. As was re- 
alized earlier I Zeldovich et al. 1968; Peebles 1968), during the recombination 
epochs photons cannot escape from the Lyman-continuum. However, at high red- 
shifts the number of neutral atoms is very small, so that the Lyman-continuum 
becomes optically thin. In order to include the es cape of photons in t he Lyman- 
continuum we follow the analytic description of Chluba & Sunyaev 12007), in 
which an approximation for escape probability in the Lyman-continuum was 
given by 



i 



1 + Tf 



(A.l) 



with rf: 



—pj j^p. Here o"i sc is the threshold photoionization cross sec- 
tion of the ls-state, Ni s is the number density of atoms in the ground state, 
and v c is the threshold frequency. For the standard cosmology the H i Lyman- 
continuum becomes optically thin above z ~ 3000 - 4000, while for He n this 
occurs at z > 12000 - 16000. As our computations show, it is crucial to include 
this process, since at high redshifts almost all loops begin or terminate in the 
Lyman-continuum (see Sect. [5}. 



A.2. High redshift solution 

At high redshifts, well before the actual recombination epoch of the considered 
atomic species, one can simplify the problem when realizing that the ionization 
degree is actually not changing significantly. Although the inclusion of intrinsic 
CMB spectral distortion does lead to some small changes in the populations 
with respect to the Saha values, the total number of electrons that are captured 
by protons and helium nuclei is tiny as compared to the total number of free 
electrons. Therefore one can neglect the evolution equation for the electrons, 
until the actual recombination epoch is entered. For H i we use this simplification 
until z ~ 3500, while for He n we follow the full system below z ~ 20 000. Before 
we simply use the RECFAST-solution for N c (Seager et al. 1999, 2000). In several 
different cases we checked that these settings do not affect the spectra. 

Furthermore, one should mention that at high redshifts for n > 2 we use the 
variable A/V, = Nj - /V 2s instead of /V;, since AiV, //V, becomes so small. Here 
N? s is the expected population of level i in Boltzmann-equilibrium relative to the 
2s-level. We then change back to the variable Nj at sufficiently low redshifts. 



A.3. Recombination and photoionization rates 

The computation of the photoionization and recombination rates fo r many levels 
is rather time-consuming. In an earlier version of our code jChluba et al. 2007a) 
we tabulated the recombination rates for all levels before the actual computation 
and used detailed balance to infer the photoionization rates. This treatment is 
possible as long as the photon and electron temperature do not depart signifi- 
cantly from each other, and when the background spectrum is given by a black- 
body. Here we now generalize this procedure in order to account for the small 
difference in the electron and photon temperature, in particular at low redshifts 
(z i 800), and to allow for non-blackbody ambient photon distributions. 



At high redshifts (z > 3000) t he electron temperature is alw ays equal to the 
Compton equilibrium temperature (Zeldovich & Sunyaev 1969): 



7/f : 



k^Ty Jx 4 n y (x)dx 
h 4 § x 3 n y (x)dx 



(A.2) 



within the given ambient radiation field. Due to the extremely high specific en- 
tropy of the Universe (there are ~ 1 .6 X 10 9 photons per baryon) this temperature 
is reached on a much shorter time-scale than the redistribution of photons via 
Compton scattering requires. For a /.(-type distortion Tp is always very close to 
the effective photon tem perature, whereas for a y-t ype distortion with y <K 1 one 
has T c » T r [l + 5.4 v] llllarionov & Syunyaev 1975a). This simplifies matters, 
since there is no need to solve the electron temperature evolution equation, and 
the photoionization and recombination rates can therefore be pre-calculated. 

At redshifts below z. ~ 3000 we solve for the electron temperature account- 
ing for the non-blackbody ambient radiation field. In this case the photon tem- 
perature inside the term due to the Compton interaction has to be replaced by 

as given by Eq. IA.21 . such that the temperature evolution equation reads 



dz H(z)[\ + 



Z] 1+/He+X, 



[Tc - ci ■ 



27/ c 
~z' 



(A.3) 



where at c = 4.91 X 10~ 22 s _I K 4 . 

Since for small intrinsic CMB spectral distortions the correction to solution 
for the temperature of the electrons is rather small, it is always possible to use the 
standard Recfast solution for T c as a reference. Tabulating both the photoioniza- 
tion and recombination rates, and their first derivatives with respect to the ratio 
of the electron to photon temperature p = T e /T r , it is possible to approximate 
the exact rates with high accuracy using first order Taylor polynomials. To save 
memory, we only consider all these rates in some range of redshifts around the 
current point in the evolution and then update them from time to time. At high 
redshifts we typically used 200 points per decade in logarithmic spacing. At 
low redshift (z < 5000) we use 2 points per Az = 1. Another improvement 
can be achieved by rescaling the reference solution for T c with the true solution 
whenever the tabulated rates are updated. With these settings we found excellent 
agreement with the full computation, but at significantly lower computational 
cost. 



A.4. Inclusion of electron scattering 



As ha s been mentioned bv l D ubrovich & Stolvarovl i 19971) and shown more de- 
tail bv lRubino-Martin et 311 12007), the broadening due to scattering of photons 
by free electrons has to be included for the computation of the He n recombina- 
tion spectrum. Similarly, one has to account for this effect, when computing the 
spectral distortions arising from higher redshifts. Here we only consider redshifts 
z < 5 X 10 4 , and hence the electron scattering Compton- y-paramete£3 



y e (z) 



=1 



:N e o-j 



m e c- 



H{z')(\ + Z') 



dz *4.8x 10~ n [1 +zf 



(A.4) 



is smaller than ~ 0.12, so that the line-broadening due to the Doppler-effect is 
significant (Ay/y|D ppi cr ~ 0.58), but still rather moderate in comparison with the 
width of the quasi-continuous spectral features, arising from high redshifts (see 
Sect. l5.ll and l5.2t . However, already at z. < 2.5 X 10 4 one has y c < 0.03, such that 
Av/v| Doppli:r < 0.29. 

Regarding the line-shifts due to the recoil-effect one finds that they are not 
very important, since even for the Hi Lyman-continuum one has Av/v| reco ii 
-0.14 at z < 5 X 10 4 . Although for the He n Lyman-continuum the shifts due to 
the recoil-effect is four times larger, we shall not include it in our results. One 
therefore expects that at frequencies v > 57 GHz the presented distortions may 
still be modified due to this process, but we will consider this problem in a future 
paper. 

For the bound -bound spectrum we follow the procedure described in 
iRubino -Martin et a^ 120071) . where the resulting spectral distortion at ob- 
serving frequency v for one particular transition is given by (see also 
Zeldovich & Sunvaev 1963) 



Mij(y)\ 



r v 3 A/, 7 (v ) 

Doppler J y 3 fl— - 



(ln[v/v ]+3v c )- 
4)i 



dvp 
vn 



— (AS) 



6 Note that y c differs from y as defined in Eq. Q}, since it describes the 
redistribution of some photon over frequency because of electron scat- 
tering rather than the global energy exchange with the ambient black- 
body radiation field. 
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Here Aljj(vo) denotes the spectral distortion for the considered transition evalu- 
ated at f requency vp and computed without the inclusion of electrons scattering 
(e.g. see Rubino-Martm et al. 2006), but accounting for the non-blackbody am- 
bient radiation field. Note that v c (z) has to be calculated starting at the emission 
redshift z cm = Vy/vo _ 1. where V;j is the transition frequency. 

For the spectral distortion resulting from the free-bound transitions one in 
addition has to include the frequency-dependence of the photoionization cross 
section. We shall neglect the line broadening because of electrons scattering for 
the moment. Then, following Chluba & Sunyaev ( 2006) and using the definitions 
of Sect. 13.21 in the optically thin limit the spectral distortion of the CMB at 
observing frequency v due to direct recombinations to level (' is given by 



A7 fc (v) 



2/iv 3 r° 



(v z ,z) 



riV,o-,(y : ) 
H(z)(l + z) 



dz, 



(A.6) 



with v- = v(l + z), 1 + Zt = Vic/v an d x z = hv z /kT y = hv/kTo. Furthermore, 
n 7 (v z , z) denotes the intrinsic CMB occupation number at redshift z including the 
spectral distortion and evaluated at frequency v z . For the Lyman-continuum one 

pLy-c, 



would in addition multiply the integrand of Eq. )A.6i by P c ^. (z), to obtain the 
approximate solution for the resulting distortion. 

Now, to include the broadening because of scattering by electrons one has 
to solve the 2-dimensional integral 



2/!V 3 

A/,cO)|Doppler = J 



where v z = v (1 + z) and 
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(A.7) 



(A.8) 



In the numerical evaluation of these integrals it is advisable to make use of the 
knowledge about the integrand, since otherwise they may converge very slowly. 



A.5. Estimate regarding the free-free process 

The free-free optical depth, Tg , is given by 

N e <x T c dz' 



r ff (x, z, z f ) 



f' K s (x,z') 



H(z')(\+z') ■ 



(A.9) 



Here N e is the free electron number density and H(z) is the Hubble factor, which 



in the radiation dominated era (z > 3300) is given by H(z) ss 2.1 X 10 
z] 2 s . The free-free absorption coefficient Kg(x, z) is given by 



'[1 + 



a At 



2n VBff x 3 ft* 



[1 - e~* 
3 (,V2 



N b gf(x,8 r ), 



(A. 10) 



where A e = h/m s c = 2.426 X 10 10 cm is the Compton wavelength of the elec- 
tron, a ss 1/137 is the fine structure constant, gjj! is the free-free Gaunt factor 
for hydrogen. We also introduced the dimensionless temperature of the photon 
field By = k^Ty/m^c 2 ss 4.6 X 10~ 10 [1 + z]. For simplicity we again shall as- 
sume 7" c = T y . Furthermore we approximated the He ++ free-free Gaunt factor 

by ?|f C++ ~ ^Sff.p an d assumed that z ^ 8000, since in the considered frequency 
range most of the free-free absorption occurs well before He in — > He II recom- 
bination (see below). 

Using the condition Tg ss 1 one can estimate the frequency xg(z), below 
which one expects free-free absorption to become important. Since at z ^ 8000 
all the atoms are ionized, the number density of free electrons is given by N e = 
(1 - Kp/2) N b ss 2.2 X 10~ 7 (1 + z) 3 cm 4 . Here we used N bi0 = 2.5 X 10~ 7 cm~ 3 
as the present day baryon number density. For xg "*C 1 one then finds Tg (x, z = 
0,z em ) ~ 9.6 x 10- 7 [g ff /5] ^T+zl^x- 1 , where 

Zem is the redshift of emission. 
Here we are only interested in photons that can be observed at x > 0.02, i.e. y > 
1 GHz today. At this frequency Tg ^ 1 for z em :> 1-7 X 10 5 . Below this redshift 
one can neglect the free-free process in the computation of the bound-bound and 
free-bound spectra. However, a more complete treatment will be presented in a 
future paper. 

Assuming that z em ~ 8000 one finds that T ff < 0.21 for x > 0.02. This 
justifies the approximations that were made above, since the contributions to the 
free-free optical depth coming from z S 8000 are not very large. 



Appendix B: Analytic solution 

B.1. The 2-shellatom 

Including only 2 shells one can analytically derive the solution for the Lyman- 
a line under quasi-stationary evolution of the populations. For this we need to 
determine the net radiative rate, Ai?L yo = A21U + M21) 7 Vap[l - wNi s /N2$ x 
n 2l/(l + "2l)L with w = 3 and H21 = «y(V2l. ^V)- ^ e rat ^° A = «2i/(l + «2i) 
is directly determined by the given ambient radiation field including the spec- 
tral distortion. Since the distortions are considered to be small we can use 
A21U + »2l) 7 V2p ~ ^2l(l + n 21^2p ^ or ^ e term m front of the brackets. 
Here and are equilibrium values for the photon occupation number and 
the 2p-population, respectively. Therefore we only have to determine the ratio 
^ = wJVis/JVzp in order to compute the Lyman-o- line intensity analytically. 

We shall first consider the situation in the case of hydrogen at high redshifts 
(z ^ 3000 - 4000). There the escape probability in the Hi Lyman-a line and the 
Hi Lyman-continuum are very close to unity. Therefore the 2s- 1 s-two-photon 
transition does not play any important role in defining the number density of 
atoms in the ground state. Furthermore, one can assume that the 2s-population is 
always in Saha-equilibrium with the continuum, and hence /V2 S ss N C N V ff2s //?2s 
where even N p ss JV H , since the total fraction of neutral atoms is tiny. 

For the Is- and 2p-states the rate equations read 

N,N„ tt U -—N 2p + A 2 i(l + n 2 l)iV2p[l - fA] * (B.la) 

w 

N,N p a 2p -/? 2p N 2p - A 2 i(l + n 2l )N 2p l\ - fA] « 0, (B.lb) 
where we substituted N[ s = t;N 2p /w and A = «2l/(l + «2l)- Solving this system 
tlfl/Vp ~/V H forf one finds 

_ a\ s P 2p +A 2 i(a ls + a 2p )[l +n 2 l] ^ ^ 

<X2 V Pls/w + Mi(an +«2p)«21 
With appropriate replacements the same expression can be used to compute the 
He 11 Lyman-o- line. 

B.1 .1 . Including the Lyman-ff and continuum escape 

In order to include the escape probability in the Lyman-ff line, P 2 \, and the 
Lyman-continuum, Pi c , one simply should replace A21 — » P 2 \A 2 \, a\ s — » 
Plcffis and ji\ s — » Pio/Sis, where the escape probabilities can be computed using 
equilibrium values for 7Vi s and N 2p . As long as the 2s- 1 s-two-photon transition 
can be neglected this yields a very accurate approximation for the Lyman-cc line 
(cf. Sects.Hm. 

Around the region where the Lyman-continuum is becoming optically thick 
(z ~ 3000 for H 1 and z ~ 1 1000 for He 11), for simple estimates one can use 

csc ^ (1.2 X 10~ 24 e v >> [1 + z] 4 for Hi 

T Ly-c ~| 3 lx lQ-24gX, s [J + z ]7/2 f or He „ 



(B.3) 



3.0 x 10 
5.0 x 10 



e v ^[l+z] 3 for Hi 



(B.4) 



with xu ~ 5.79 x 10 4 Z 2 [l + zi~ 



B.1.2. More approximate behavior 

In order to understand the solution for the H 1 Lyma n-g fine we now turn to 
the corresponding intensity as a function of redshift I Rubino-Martm et ^120061 
e.g.). This yields 
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An H(z)[\ + z] 3 An H(z)[l + z] 3 

Using the approximation (B.2\ for £ one can then find 

a2pPlJw(\ + n 2 \) - a ls /3 2p n 2i 



l-fA< 



[Q-2p/?ls/vf +A2l(ffl s +Q-2p)H2l](l +«2l) 



(B.5) 



(B.6) 



With this one then has 



he A 2l N 2p 

My S3 — 



Q-2p/3ls(l +«2l) 



An H(z)[\ + z] 3 a 2p /3i s + wA 2 i(Q!i s + Q-2p)i2i 



-fel+Kl) 



(B.7) 



Here we used ^21 = hv 2 \jkTy and (1 + «2l)/«2i = e X2,+ni , with the frequency 
dependent chemical potential /J21 = fi(X2l)< 



1 Note that even if one (more correctly) uses N p = Nn-Ni s -N 2 s-N 2p 
in Eq. dB . 1 1 > the solution for £ does not change. 
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First factor 

We can now simplify the expression )B.7t when realizing that except for 
the term inside brackets, in the case of small intrinsic CMB spectral distor- 
tions, one can just use equilibrium values. At high redshifts one has H(z) cc 
(1 + z) 2 . Furthermore 1$ * 3N^e~ x u » 3 N e N H a\ q s <r* 21 /ff* and 0: q = 

orT* e - "'" I fi(T e ). Also with rather high accuracy one finds as and 
^Epi <g A 2 \(ai s + a 2p )n 2 i, so that 

_ he A 2I N 2p a2 P /Si s (1 + n 2 l) 



to H(z)[l + z] 3 a 2p /3n + wA 2 i(a u + a 2f )n 2 i 

hc_ 3N e N H a 2p (1 +z) m 
to ff(z)[l + z] 3 4 + r[- c K 4 + rf 1 c ' 



(B.8) 



Here we have also included the escape probabilities in the Lyman-a line and 
continuum as explained in Appendix |B.l.l| Note that the Lyman-a escape prob- 
ability drops out of the expression, so that only the Lyman-continuum escape 
probability is strongly affecting the pre-recombinational line shape. We have also 



used a? = ^fi(T c ) e x " /' q , with the integral 

? = £ dv * <Ti(Vfc) v - M -' (Xic) ' 



(B.9a) 



where we have assumed cr,(v) ss cr,-(Vj c ) Jf-. The integral M{(x) is defined and 

discussed in Appendix [C] For the 2p-state one has o"2p(v2pc) v^, ~ 7.54 X 

10 27 Z 4 cm- 2 s- 3 . 

We checked the scaling of F numerically and found 



F(z) * 5.6 x 10- 26 , (1 \f 2 IA Jm- 2 s- 1 Hz" 1 sr" 1 

1 +T Ly-c/ 4 

within < 20% accuracy in the important redshift range. 



(B.10) 



Second factor 

Using the definitions of a, and p 1 , as given in Sect. l3.1l and l3~2l one directly finds 



w e-^i+wi) = e -m G(z) with 



G(Z) : 



) ls <">2p 



Here Ap = 1 - p and we introduced the notation 



</(y)>, = vV,(v)/(v)dv 



for the average of some function /(v) over the photoionization cross-section of 
level i. 

In full thermodynamic equilibrium one has G ctl (z) = 1, a property that 

can be verified using Eq. fBTfl with n = and p = 1, since then 

| Be l*)+(%-*/py = (« pl ).. Therefore we can write G = 1 + AG. Using 
(/>; = < / ec i >, + < A/ >,, for small intrinsic CMB distortions (i.e. < Af >,/< / >, « 
1 ) one finds 
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2 P 



("Pl)ls ("Pl) 2 p 
where ff! = (x- A", c )Ap/p. Putting things together we then have 
i 

l-w ■ 



(B.13) 



«2p/3ls 



("P'^p ("P') ls 

(»pi^ P ) 2p (»piO, 



("Pl) 2p ("Pl)ls 



i p 2 l +P2pc -Plsc , 



(B.14a) 
(B.14b) 



To lowest order, Eq. )B.14b> shows that the main reason for the emission in 
the Lyman-ff line is the deviation of the effective chemical potential from zero 
at the Lyman-a resonance, and the Lyman- and Balmer-continuum frequency. 
However, the averages over the photoionization cross-section still lead to some 
notable corrections, so that also the small difference in the electron and photon 
temperature plays a role. 

If we again use the Kramers-approximation for the photoionization cross- 
v? ._. 
section, o",(v) » o"((v; c )-y, looking at Eq. Qj for p in the case of a small y-type 

distortion, one can write 



(dpi ). ss KjM-i(x ic ) 
( « p i p ). « Kj y [4 MoO; c ) - 5 (x, c )] 

(«pl^).*'- V| 



-[M ( 



(B.15a) 
(B.15b) 

(B.15c) 



where ff, = const and the integrals S and Mq are defined in Appendix|C] Keeping 
only the leading order terms, we have 



l-w 



aup 2p 
«2p/Sls " 

-yxisc [6.3 -0.9375 ,v lsc +4.7 e 
>.4 - x 2 



1.175 <r*i«/ 4 ] 



+ 1 



(B.16a) 



It is important to mention that this is still a rather rough approximation, since 
already applying the Kramers-formula for the photoionization cross-section in- 
troduces some significant simplification. However, this approximation may be 
useful for simple estimates. 



Appendix C: Some integrals 

C.1. Integrals M/ 

In the evaluation of the recombination and photoionization rates, integrals of the 
form M, = f x 1 dx/[e x — 1] appear. Below we now discuss those which are of 
importance for us here. 



C.1.1. Integral M. 



(B.ll) Fori: 



(B.12) 



- 1 one can write 
dx 



£ 



y,Ei(te ic ) k 



hv ic <kT y 
I 1 



11 -6y 
12 " 



J2 

ln(A, c ) 

12 2 v ' 



(C.la) 



(C.lb) 



where y ss 0.5772 is the Euler constant and we made used of the exponential 
integral Ei(jc) = J e~' dt/t. In the limit hv\ c < kT y the given approximation is 
accurate to better than 1%. For hvj C > kT y the first five terms in the full sum also 
yield similar accuracy. 

Since x ic a 5.79 X 10 4 Z 2 «- 2 [1 + z]" 1 , it is clear that at z < 1.45 X 10 4 Z 2 
both Lyman- and Balmer-continuum are still in the exponential tail of the 
CMB blackbody. In the redshift range 1.45 X 10 4 Z 2 < z < 5.79 X 10 4 Z 2 
the Lyman-continuum is still in the exponential tail of the CMB, while the 
Balmer-continuum is already in the Rayleigh- Jeans part of the spectrum. Only at 
z > 5.79 X 10 4 Z 2 one can use the low frequency expansion of Eq. jB.9t for both 
cases. However, to within < 30% one may also apply Eq. IC.lal in the whole 
range. 



C.1.2. Integral M 

For i = one can write 



M 



: f" dx/[e x - 1] = £ r e* dx = f j e- kx - 



Ik 



-x, c -ln(e x '< -1), 



(C.2) 



which for x/ c < 1 can be approximated as Mq ~ %r ~ 327 - lnfec), while for 
Xic » 1 one has Mq « e _rit [1 + e~*fc/2]. 
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C.2. Integrals 

In the evaluation of the recombination and photoionization rates one also en- 
counters S(x) = f°° dxx e x + .U . The first part of this integral, oc xe x /[e x - l] 2 , 

can be directly taken yielding f°° dx xe x /[e x -l] 2 = jr, c £ A '> c /[<r t <' -l]-ln(e* - 1). 

Jx ic 

Introducing the polylogarifhm Li„(;t) = x^/k" and realizing x/[e x - l] 2 = 
Y,k =l kxe- {k+l)x one can find 

e' c — 1 

- (2 - x ic )]n(e x * - l)-Li 2 (<T* fc ) 

hv jc >kT Y mK g 

" Znr 11 *'* 1 ^ (C3a) 



